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摘要 :准确 地 划分 地 带 性 主要 植被 类 型 分 布 的 适宜 区 域 , 可 为 区 域 植被 恢复 与 重建 .生物 多 样 性 保护 等 工作 提供 有 着 的 理论 参 
考 。 在 检验 1961 一 2013 年 东北 地 区 气候 变化 突变 点 的 基础 上 ,基于 东北 地 区 主要 植被 类 型 热量 指标 ,人 研究 气候 变化 对 该 区 域 
主要 植被 类 型 适宜 分 布 区 域 的 影响 ,并 利用 2000 年 和 2013 年 土地 类 型 覆盖 数据 (MCD12Q1) ,对 推算 的 适 攻 分布 区 域 进行 验 
证 ,结果 表明 :大 兴安 岭 气温 突变 点 为 1982 年 ,其它 地 区 为 1988 年 ,东北 地 区 各 生态 地 理 区 年 降水 量 没有 明显 的 突变 点 。 气 温 
突变 点 后 适宜 在 东北 地 区 生长 的 主要 植被 种 类 没有 变化 ,但 各 植被 类 型 的 分 布 区 域 均 有 所 改变 ,其 申 高 山 冻 原 、 亚 高 山 矮 曲 林 、 
寒 温带 针 叶 林 和 温带 针 阔 叶 混 交 林 适宜 分 布 区 域 面积 减 小 , 暧 温带 落叶 阔 叶 林 和 温带 莘 原 面积 增 大 ,Ss 气温 突变 点 前 后 各 植被 
类 型 适宜 分 布 区 的 地 理 中 心 均 发 生 了 不 同 程 度 的 移动 ,其 中 移动 距离 最 大 的 是 南部 地 区 的 亚 高 山 矮 曲 林 ,向 东北 方向 移动 了 
135.44km。 与 主要 植被 类 型 实际 分 布 相对 比 ,证 明 气候 变化 对 研究 区 植被 类 型 分 布 名 能 已 经 产生 了 影 
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Abstract: A suitable distribution of zonal vegetation types could provide a reference for the restoration and reconstruction of 
regional vegetation and protection of biodiversitys Between 1961 and 2013, we used mutation point detection and the thermal 
index to study the impacts of climaté \change on the distribution of the main vegetation types in Northeast China. The land 
cover data for 2000 and 2013 fromNASA LPDAAC (The Land Processes Distributed Active Archive Center) EOS were 
used to verify the simulated distiibution regions. The results showed that the temperature mutation point for Daxing' an 
Mountains was in 1982% anduthe mutation points for the other regions were in 1988. The annual precipitation at each eco- 
geographical region@showed no obvious mutation. After 1982, the main vegetation types, which were suitable for growth ， 
showed no obvigis changes’; however, their distribution regions had been changed. The areas of suitable distribution regions 
for alpine tundra, subalpine forest, cold-temperate coniferous forest, and temperate mixed forest decreased continuously; 
however, the areag for warm-temperate deciduous broad-leaved forest and temperate grassland increased. The geographic 
centers for the suitable distribution of each vegetation type had moved to different degrees before and after the temperature 
mutation points. Among them, the center of the subalpine forest had moved towards the northeast, and the distance was the 
longest (135.44km). On the basis of the comparative analysis, we concluded that climate change possibly had an impact on 


the distribution of vegetation types in Northeast China. 
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全 球 气 候 变 化 及 其 带 来 的 影响 已 受到 广泛 关注 ,IPCC 第 五 次 评估 报告 指出 ,2003 一 2012 年 平均 温度 比 
1850 一 1900 年 上 升 了 0.78% ,而 降水 在 地 区 和 季节 上 的 差异 也 有 所 增加 ,未 来 极端 性 天 气 气候 事件 的 发 生 概 
率 可 能 将 进一步 增加 所 。 中 国 的 气候 变化 与 全 球 的 变化 趋势 是 一 致 的 ,平均 增 温 速 度 高 于 全 球 或 北半球 同 
期 ,而 北方 和 青藏 高 原 尤 为 明显 '* ;降水 趋势 变化 存在 明显 的 区 域 差 异 * 引 ,其 中 东北 地 区 ' 和 西部 地 区 是 
的 降水 出 现下 降 趋势 。 

气候 作为 生态 系统 演变 的 重要 控制 因子 , 它 的 任何 变化 都 会 对 生态 系统 产生 影响 ”中 ,进而 对 地 带 性 植 
被 的 建 群 种 分 布 产生 影响 "中 ,所谓 地 带 性 植被 是 指 在 分 布 上 与 气候 带 界 线 大 体 相符 的 大 面积 植被 类 型 , 它 
综合 反映 了 各 自 的 大 气候 特点 和 一 定 的 生态 地 理 空间 王国 。 有 学 者 沿用 Kira 的 温暖 指数 (WI 和 赛 冷 指数 
(CIT) 划分 地 带 性 植被 !55] ,有 学 者 以 WI 为 基础 进行 拓展 ,或 在 此 基础 上 引入 湿度 指数 (HID) ,这 些 划分 
方法 均 是 建立 在 气候 是 地 带 性 植被 类 型 分 布 的 主要 影响 因子 基础 上 的 ,因此 对 东北 地区 的 划分 基本 一 致 , 即 
东北 地 区 的 自然 植被 属 欧 亚 森林 草原 植物 亚 区 和 中 国 - 日 本 森林 植物 亚 区 , 蔬 要 地 带 性 植被 为 寒 温 带 针 叶 
林 、 温 带 针 叶 混 交 林 和 暧 温带 落叶 阔 叶 林 。 东 北 区 地 处 全 球 气 候 变 化 的 敏感 区 域 , 地 上 开 辽 阔 , 地 势 复 杂 , 气 候 
变 率 大 9, 因此 许多 学 者 就 气候 变化 对 该 区 域 植被 分 布 产生 的 影响 做 了 一 些 研究 二 如 冷 文 芳 等 .在 模拟 
100a 后 气温 和 降水 变化 的 基础 上 ,使 用 logistic 回归 模型 预测 东北 森林 主要 建 群 树种 的 空间 分 布 变 化 ,结果 表 
明 多 数 树种 的 覆盖 率 可 能 降低 ,只 有 红 松 .长白 落叶 松 和 蒙古 栎 兴 可 能 增加 。 程 肖 侠 ' 站 也 在 预测 模拟 未 来 
100a 气候 情景 下 东北 地 区 森林 类 型 的 变化 ,认为 气候 变 暧 不 利于 东 毕 主要 森林 类 型 的 生长 ,主要 的 针叶树 种 
比例 可 能 会 减少 ,阔叶树 种 比例 会 增加 。 以 上 研究 均 是 建立 在 气候 预测 基础 上 的 ,但 目前 研究 表明 东北 地 区 
的 气候 变化 已 成 不 争 的 事实 ,这 种 变化 对 该 区 域 森林 主要 建 群 树种 的 分 布 区 域 是 否 已 经 产生 了 影响 ,影响 
程度 如 何 还 少见 报道 ,本 研究 正 是 针对 这 一 问题 往 总 结 前 人 研究 成 果 的 基础 上 ,收集 近 十 几 年 新 增加 的 科研 
资料 ,根据 东北 地 带 性 建 群 树种 和 常见 种 温暖 指数 范围 ,从 生态 气候 学 的 角度 探讨 气候 变化 对 东北 三 省 森林 
植被 地 理 分 布 的 影响 。 该 研究 不 仅 有 利于 民 域 森林 植被 的 恢复 与 重建 ,也 可 为 该 区 域 自然 资源 的 合理 利用 、 
土地 退化 防治 .生物 多 样 性 的 保护 等 区 域 可 持续 发 展 战略 和 规划 的 制定 提供 理论 参考 。 


1 研究 区 概况 


东北 三 省 包括 黑龙 江 省 吉林 省 和 辽宁 省 ,位 于 118.83° 一 135.09°E,53.56° 一 38.72°N 之 间 , 南 邻 黄 海 和 
渤海 , 东 面 和 北面 有 蝎 绿 江 图 们 江 . 乌 苏 里 江 和 黑龙 江 环 绕 , 西 面 与 内 蒙古 接壤 ,总 土地 面积 78.7 万 km?, 境 
内 有 大 小兴安岭 长 向 山系 和 东北 平原 (包括 松 辽 平原 .辽河 平原 和 三 江平 原 ) ,海拔 高 度 在 2 一 2667m 之 间 。 
该 地 区 属 温带 季风 气候, 其 中 北部 大 兴安 地 区 属 寒 温 带 气候 ,其 它 地 区 属 温 带 气候 , 年 平均 气温 -4.15 一 
10.99%C ,年 降水 量 366.90 一 1080.47mm ,日 平均 日 照 时 数 5.22 一 8.10h。 


2 数据 来 源 与 处 理 


2.1 号 理 信息 数据 

数字 高 程 模 型 (DEM) 为 SRTM 地 形 产品 V4.1 版 本 数据 ,空间 分 辩 率 为 90m ,数据 来 源 于 中 国 科 学 院 计算 
机 网 络 信息 中 心 国际 科学 数据 镜像 网 站 (http://www.gscloud.cn) 。 

研究 所 需 的 县 级 和 省 级 行政 区 划 数 据 来 自 于 中 国 气 象 局 下 发 的 1:25 万 基础 地 理 信息 ,对 数据 进行 拓 补 

查 , 去 除 省 界 和 县 界 的 间 际 。 地 面 观测 站 位 置 数据 为 中 国 气象 局 下 发 的 矢量 数据 。 
2.2 土地 覆盖 类 型 数据 

土地 类 型 覆盖 数据 (MCD12Q1) 来 源 于 美国 NASA LPDAAC(The Land Processes Distributed Active Archive 
Center) EOS 数据 中 心 ,该 数据 的 空间 分 辩 率 为 S00m ,是 利用 当年 Terra 和 Aqua 星 观测 所 得 的 数据 描述 十 地 
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窗 盖 类 型 ,信息 提取 主要 技术 是 监督 决策 树 分 类 ,产品 包含 5 种 不 同 土地 覆盖 分 类 方案 的 分 类 结果 和 对 应 的 
质量 信息 数据 集 , 本 研究 选用 “IGBP 的 全 球 植被 分 类 方案 ”的 分 类 结果 ,该 分 类 方案 根据 国际 地 圈 生 物 圈 计 
划 (IGBP) ,本 文 所 用 数据 的 时 间 为 2001 年 和 2013 年 ,空间 范围 为 东北 地 区 ,利用 MODIS Reprojection Tool 
software( MRT) V4.0, 进 行 影像 的 拼接 、 重 采样 (统一 空间 分 辨 率 为 500m) 和 投影 转换 (由 Sinusoidal 投影 转换 
为 经 纬度 投影 ) 等 预 处 理 。 
2.3 气象 数据 

气象 数据 来 源 于 黑龙 江 省 .吉林 省 和 辽宁 省 气象 局 
的 1961 一 2013 年 的 196 个 气象 台 站 逐日 整编 资料 , 包 
括 日 平均 气温 和 日 降水 量 , 站 点 的 地 理 位 置 和 海拔 高 度 
如 图 1。 
2.3.1 气象 数据 预 处 理 

气候 变化 研究 必须 建立 在 可 靠 的 资料 基础 上 ,而 由 
于 气象 台 站 迁移 .观测 数据 误差 和 缺失 城市 热岛 效应 
等 原因 ,会 产生 数据 序列 的 不 均一 性 ,因此 在 进行 分 析 > 
之 前 需要 对 数据 进行 预 处 理 , 最 大 限度 地 降低 因数 据 误 。 4 ” ee 
差 而 影响 分 析 结 果 的 准确 性 。 E 

(1) 质 量 控制 120° 125° 130° 135°E 

由 于 观测 仪器 、 资 料 处 理 和 传输 过 程 可 能 引起 误 
差 ,资料 中 很 难 避免 存在 一 定 的 错误 ,最 大 限度 地 减少 Fig;l_Distribution diagram of meteorological observation 
这 些 错 误 的 可 能 影响 是 非常 必要 的 ' 引 ,因此 本 研究 采 /sthtions and eco-geographical in the Northeastern China 
用 翟 盘 茂 ! 约 的 方法 对 气象 数据 进行 质量 控制 ,删除 错 “和 :大 兴安 岭 . 工 : 松 辽 平原 中 部 、 开 :东北 东部 山 前 平原 .TV : 东 
误 数 据 。 另 外 Karl525] 也 指出 ,把 与 时 间 有 关 的 资料 缺 北 东部 山地 、V ;三江 平原、V[: 华 北山 地 丘陵 、 讶 :辽东 股东 山地 
测 假设 为 零 时 ,会 在 气候 变化 趋势 研究 中 丙 起 误 丑 本 。 区 
研究 中 监测 缺 测 数据 的 方法 是 当 一 年 内 缺 测 超过 20d 
( 约 大 于 一 年 内 实际 观测 日 数 的 59%) 时 ,该 年 的 统计 量 被 设 为 缺 测 , 当 45a 序列 中 缺 测 总 量 超 过 5a 时 , 则 该 
测 站 被 排除 使 用 。 

(2) 正 态 分 布 检验 

多 数 气候 诊断 方法 和 预测 模型 是 在 气候 变量 呈正 态 分 布 假定 前 提 下 进行 的 ,所 以 对 气候 变量 是 否 呈 正 
分 布 形 态 的 检验 是 十 分 必要 的 。 本 文 利用 SPSS 软件 对 各 气候 区 逐日 气温 和 相对 湿度 数据 进行 正 态 检验 ,4 
著 性 水 平 sig. 有 98.7% 均 三 0.05 ,说 明 各 气候 区 逐日 气温 数据 绝 大 多 数 近似 正 态 分 布 ,分 析 时 可 不 对 其 进行 标 
准 化 处 理 。6 与 气温 数据 不 同 , 日 降水 量 的 变化 并 不 具有 逐步 的 连续 的 特征 ,因而 判断 降水 数据 是 否 服从 正 态 
分 布 是 以 月 为 自 变 量 进 行 检 验 , 检 验 结果 绝 大 多 数 服从 正 态 分 布 , 因 而 分 析 时 不 作 标 准 化 处 理 。 

(3) 热 锅 效 应 分析 

IPCGs 3 次 评估 报告 认为 ,从 全 球 平均 来 看 ,城市 化 对 地 表 平 均 气 温 记 录 的 影响 微乎其微 ,其 实际 作用 
要 比 观测 到 的 增 温 值 小 一 个 数量 级 '*” ,而 Wang'” 和 赵 宗 蓄 ' 池 却 认为 ,热岛 效应 对 长 期 气温 序列 研究 的 影响 
不 容 忽 视 , 于 是 本 研究 根据 2015 年 统计 年 鉴 中 人 口 数据 ,删除 人 口 大 于 等 于 100 万 人 和 人 口 大 于 等 于 50 万 
人 的 城市 测 站 的 气温 观测 数据 ,计算 未 删除 .删除 人 口 大 于 等 于 100 万 人 城市 和 删除 人 口 大 于 等 于 50 万 人 的 
城市 的 年 平均 气温 变化 速率 ,结果 分 别 为 0.0149 .0.0126 .0.0145% va ,它们 之 间 的 差别 比 增 温 值 (21 世纪 以 来 
比 19 世纪 60 年 代 高 0.16% ) 小 一 个 数量 级 ,因此 本 研究 不 考虑 城市 热岛 对 气温 序列 变化 的 影响 。 

(4) 数 据 插 补 

本 文 对 缺 测 测 站 的 气象 数据 采用 Newton 法 进行 搬 补 。 


图 1， 东 北三 省 气象 观测 站 和 生态 地 理 区 分 布 示意 图 
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2.3.2 ”气象 数据 突变 检测 


通过 对 东北 地 区 1961 一 2013 年 气温 和 降水 突变 点 的 检验 ,为 探寻 和 气候 变化 对 研究 区 主要 地 带 性 植被 类 


型 分 布 可 能 带 来 的 影响 提供 时 间 节 点 。 
(1) Mann-Kendall 检验 
Mann-Kendall( 简称 MK) 是 气象 学 /气候 学 中 经 常用 来 进行 突变 检验 的 一 种 方法 ,具体 方法 如 下 : 
设 原始 时 间 序 列 为 y, ,y,,…,y, ,mi 表示 第 i 个 样本 y, 大 于 y(1 三) 短裤 的 累积 数 ,定义 统计 量 : 


k 
di= > m(2<ksn) 
t= 


在 原 序列 随机 独立 等 假设 下 ,d 的 均值 和 方差 分 别 为 ， 
k(k—1) 
E(d) = 一 
rd) hs- D (2k+ 35) 
ee 72 
将 上 面 公式 的 标准 化 ,得 ， 
di 一 E(d;) 
UF, = 
/varl d,) 


UF 组 成 一 条 UF 曲线 ,通过 信 度 检验 可 得 出 其 是 否 有 明显 的 变化 趋势 。 把 此 方法 引用 到 反 序 列 中 ,计算 


条 
得 到 另 一 条 曲线 UB, 则 两 条 曲线 在 置信 区 间 内 的 交点 确定 为 突变 点 。 


a 才 \ 
滑动 1 检验 法 , 它 是 用 来 检验 两 随机 样本 平均 值 的 显著 性 差异 ”1 。 为 此 ,把 一 个 长 度 为 n 的 连续 随机 变 
量 % 分 成 两 个 样本 子 集 x 和 %%,, 让 jw、5s 和 nn 分 别 代表 名 的 平均 值 样本 方差 和 样本 长 度 (i=1,2)。 两 样本 子 


(2) 滑 动 了 检验 
集 始终 间隔 一 个 样本 ,这 样 检 验 的 就 是 某 一 年 后 太 年 和 前 n, 年 均值 的 显著 性 差异 。 


原 假 设 已 :LI 一 =0。 定义 1 统计 量 为 : 


(CO -1)S7 +( — 1)5? 


这 里 ,5, 是 联合 样 术 方差 ， 
S, ni+n,—2 
为 o, 的 无 偏 估计 (E[ Sp ] =o” ) ,显然 如 -i(ni+ns-2) 分 布 ,给 定 信 和 度 a, 得 到 临界 值 i, ,计算 6 后 在 H, 下 
著 性 差异 。 当 1i61<i, 时 , 则 接受 原 假 设 肛 ,。 


EF 


MW 


比较 i。 沁 s。, we 全 .时 ,否定 原 假设 世 ,, 即 说 明 其 存在 
局 MTT 哺 检验 某 一 个 时 间 序 列 , 可 以 得 到 相应 的 统计 量 序列 ,峰值 和 谷 值 所 对 应 的 年 份 称 之 为 突变 年 


份 .突变 点 或 者 过 渡 年 份 。 
(3) 累积 距 平 检验 
累积 距 平 也 是 一 种 常用 的 、 由 曲线 判断 变化 趋势 的 方法 ,对 于 序列 x, 其 某 一 时 刻 上 的 累积 距 平 表示 为 : 


% = (x; — x) (t= 1,2,3,.…,n) 
=1 


其 中 x = 
用 
姓 


| 2 
将 nn 个 时 刻 的 累积 距 平 值 全 部 算出 , 即 可 绘 出 累积 距 平 曲线 进行 趋势 分 析 。 
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3 结果 与 分 析 


3.1 气象 数据 突变 检验 
3.1.1 研究 区 划分 

东北 三 省 跨越 寒 温 带 .中 温带 和 暖 温带 , 且 地 形 复杂 ,山地 、 丘 陵 和 平原 相间 ,海拔 高 度 跨 越 266Sm ,气候 
变化 趋势 和 突变 点 可 能 会 存在 差异 ,因此 本 研究 参照 (中 国生 态 区 划 研 究 》 2 提供 的 中 国生 态 地 理 分 区 数据 
(来 源 于 中 国 科 学 院 资源 环境 科学 数据 中 心 (http ://www.resdc.cn ) 提供 的 矢量 数据 ) ,分 别 检 验 东 北三 省 的 7 
个 生态 地 理 分 区 (图 1) 的 年 平均 气温 和 年 降水 量 的 突变 点 。 

该 生态 地 理 分 区 数据 是 基于 温度 指标 把 全 国 划 分 为 寒 温带 、 中 温带 、 暧 温带 、 北 亚热带 等 LB 个 区 域 ,根据 
干 湿 状 况 指 标 把 全 国 划 分 为 潮湿 湿润. 半 湿 润 . 半 干旱 等 6 个 区 域 ,经 上 述 2 种 指标 组 合 后 把 全 国志 分 成 46 
个 生态 地 理 分 区 ,东北 三 省 占 其 中 的 7 个 生态 地 理 分 区 (图 1)。 

3.1.2 ”突变 点 检测 

从 年 平均 气温 MK 检验 结果 可 见 (图 2) ,东北 三 省 7 个 生态 地 理 分 区 的 年 平 列 气温 均 存 在 明显 的 上 升 趋 
势 , 其 中 生态 地 理 区 工 的 UF 和 UB 在 1982 年 处 交叉 ,表明 从 1982 年 开始 该 区 的 气温 可 能 发 生 突变 ,同时 滑 
动 了 检验 结果 和 通过 计算 气温 序列 的 累积 距 均 显示 1982 年 为 该 区 后 温 突 变 年 份 ;用 同样 的 方法 检测 生态 地 
理 区 一 证 的 气温 突变 年 为 1988 年 ,该 结论 与 31 一 致 。 

年 降水 量 MK 检验 结果 与 年 平均 气温 不 同 ,生态 地 理 区 卫 、 古 VV 和 VI 的 MK 降水 量 统计 值 UF 在 [ -1.96， 
1.96] 波 动 , 说 明 这 几 个 生态 地 理 区 在 研究 时 间 段 内 年 降水 量 欧 显著 变化 ;生态 地 理 区 IV 从 1976 年 开始 显著 
下 降 ,1983 年 开始 显著 上 升 ;生态 地 理 区 三 从 1963 年 开始 下 降 , 在 2003 年 出 现 显著 下 降 趋 势 。 结 合 MK 和 滑 
' 动 了 检验 结果 ,东北 地 区 的 7 个 生态 地 理 区 年 降水 量 没有 明显 的 突变 点 ,因此 本 文 只 人 研究 气温 突变 对 东北 地 
加 区 地 带 性 植被 分 布 适宜 区 的 影响 。 

3.2 ”气候 变化 对 东北 地 带 性 植被 分 布 的 影响 
3.2.1 植被 分 布 热量 指标 WW 
按照 徐 文 铎 [的 方法 计算 水 热 指标 


温暖 指数 ( WW) : 


-5 。 式 中 ,1 为 5C 以 上 的 月 平均 气温 。 
再 根据 树种 的 地 理 分 布 资料 , 按 半 峰 宽 (PWH) 计 算法 ' 引 确定 每 个 树种 的 热量 分 布 范围 . 


-oPWH 下 限 


园 


最 适 范围 ; 局 本 
Xt PWH 上 限 
PWH=2.3545 

式 中 , 闻 为 找 量 指标 的 平均 值 ,5 为 标准 差 。 这 种 确定 方法 能 使 树种 约 有 78% 分 布 在 最 适 热 量 范围 之 内 。 
表 1 为 6 个 东北 地 带 性 植被 建 群 种 和 常见 种 植被 类 型 的 央 范围。 
3.252, “WI 插值 
把 第 3.1.2 节气 候 突 变 点 检测 结果 作为 分 界 点 ,分 别 计算 突变 点 前 、 后 各 个 气象 观测 站 点 的 忠平 均值， 
即 利 用 1961 一 1981 年 生态 地 理 区 I 和 1961 一 1987 年 生态 地 理 区 本 一 W 的 数据 计算 突变 点 前 每 个 站 点 中 的 
年 平均 值 ( WI ) ;突变 点 后 的 由 的 平均 值 ( 吧 . ) 所 用 数据 的 时 间 范 围 为 1982 一 2013 年 生态 地 理 区 工 和 
1988 一 2013 年 生态 地 理 区 卫 一 证 的 数据 ,计算 方法 同上 。 

由 于 东北 三 省 的 气象 观测 站 点 少 且 分 布 不 均 , 尤 其 在 地 形 复杂 的 山区 气象 观测 站 分 布 尤其 稀 跑 ,这 些 有 
限 的 气象 观测 站 无 法 代表 相距 较 远 地 区 的 气候 状况 ,基于 此 数据 计算 的 WW 的 代表 性 也 会 较 差 ,为 解决 这 一 
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-4 -6 -4 
1961 1977 1993 2009 1971 1979 1987 1995 2003 1961 1977 1993 2009 
年 份 Year 


2 东北 三 省 1961 一 2014 年 平均 气温 、 年 降水 量 Mann-kendall 曲线 和 滑动 了 检验 曲线 
Fig.2 Statistics curyes of Mann-kendall test and moving T-test of mean temperature and annual precipitation between 1961 and 2014 in 
Northeast China 
MK-T1 一 MKAT7 :第 工区 车 态 地 理 分 区 气温 的 MK 曲线 图 ,TTI 一 T-T7 :第 工 一 丰 生 态 地理 分 区 气温 的 滑动 了 检验 曲线 图 ,MK-P1 一 MK- 
P7: 第 [过 如 生态 地 理 分 区 降水 的 MK 曲线 图 ,T-P1 一 TP7: 第 工 一 碍 生态 地 理 分 区 降水 的 滑动 了 检验 曲线 图 ; MK 检验 和 滑动 了 检验 的 a 
=005 ,其 中 滑动 了 检验 的 步 长 为 10a 


问题 本 研究 对 单 点 的 取 数据 进行 插值 ,由 于 天 是 基于 气温 计算 而 来 ,而 气温 与 经 度 、 纬 度 和 海拔 高 度 又 具 
有 较 好 的 相关 性 '*1 ,因此 本 研究 以 经 度 、 纬 度 和 海拔 高 度 为 自 变 量 ,分 别 以 下 ,和 WW 为 因 变量 拟 合 多 元 线 
性 回归 模型 , 拟 合 结果 见 表 2。 

利用 回归 方程 推算 出 各 个 观测 站 的 ,和 由 .的 拟 合 值 ,用 拟 合 值 减 去 实际 计算 值得 到 的 差 值 分 别 记 为 
WI, 和 WL, ,再 利用 IDW 插值 方法 把 这 些 点 的 差 值 转换 为 面 ,用 这 个 面 数据 作为 订正 系数 来 订正 利用 回归 方 
程 拟 合 出 的 WW 和 WL 面 数 据 , 得 到 订正 后 的 WW 和 WI, 的 面 数据 ,再 从 中 提取 各 个 观测 点 订正 后 的 点 数据 ,并 
减 去 实际 计算 值得 到 的 差 值 分 别 记 为 天 ,和 Wi ,分 别 统计 对 比 WL 与 天 ,和 WI 与 WL 的 标准 误差 .标准 
差 和 数值 范围 ,结果 表明 , WI, 标准 误差 标准 差 和 数值 范围 分 别 比 WI 小 0.13 .1.77 .6.82cm, Wl 标准 误差 、 
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标准 差 和 数值 范围 分 别 比 WL 小 0.12 .1.57 .7.83cm ,说 明 经 过 订正 后 的 中 值 与 实际 计算 值 更 为 接近 。 


表 1 东北 地 带 性 植被 的 温暖 指数 范围 1041 


Table 1 The Scope of warmth index for the zonal vegetation in Northeast China[2] 

人 7 最 适 范围 

植被 类 型 建 群 种 和 常见 种 a 
i Rdificat J . The optimal range 

vegetation types dificators and companions of WIA(C .月 ) 
高 山 冻 原 牛皮 杜 鹏 ( Rhododendron aureum ) .仙女 木 ( Dryasoctopetala ) 、 7 15 
Alpine tundra 松 毛 染 ( Phyllodocecaerulea ) 、 圆 叶 柳 (Salix crenataHao) 
亚 豆 山 鱼 曲 林 优 松 ( Pinuspumila) 后 桦 ( Betulaermanii) 15—35 
Subalpine forest 
寒 温 带 针 叶 林 鱼 鲜 云 杉 ( Piceajezoen ) 、 红 皮 云 杉 ( PiceakoraiensisNakai)、 奥 松 »5__65 
Cold-temperate coniferous forest ( Abiessibirica) 兴安 落 叶 松 (Larixgmelinii) 、 樟 子 松 ( Pinussylvestris ) 


红 松 ( Pinuskoraiensis) 、 沙 松 ( Abiesholophylla) 、 
温带 针 阔 叶 混 交 林 赤 柏 松 ( Taxuscuspidata) 、 鹅 耳 杨 ( Carpinusturczaninowii) 、 


Temperate mixed forest 蒙古 栎 ( Quercusmongolica ) 、 紫 概 ( Tiliaamurensis) 、 39 
黑 桦 ( Betuladavurica) 
暧 温带 落叶 阔 叶 林 辽东 栎 ( Quercuswutaishansea ) 、 榭 标 〈( Quercusaliena 入 栓 皮 栎 55__95 
Tarm-temperate deciduous broad-leaved forest ( Quercusvariabilis ) 、 赤松 ( Pinusdensiflora ) 、 油 松 ( Pinustabulaéformis ) 
温带 草原 大 针 茅 ( Stipagrandis ) 、 内 加 尔 针 茅 ( StipaBaicdlensis ) 、 于 氏 针 苞 55__85 
Temperate grassland ( Stipakrylovii) . 羊 曹 (Teymauschinenszs) 、 线 叶 菊 ( Filifoliumsibiricum) 
表 2 温暖 指数 回归 方程 
Table 2 Regression equation of warmth index 
可 归 方程 举 本 炎 
WW 回归 方程 . R Sig. RMSE 样本 数 
Regression equation Sample number 
WI Y= 346.804-1.313 log-2.346 lat-0.034 alt 0.961 0.000 3.24 209 
7 7Y= 353.006-1.350 log-2.248 lat-0.037 alt 0.951 0.000 3.72 209 


7: WW,log: 经 度 (°) ,lat: 纬度 (。) ,alt: 海拔 高 度 (m) 


3.2.3 ”植被 类 型 分 布 区 域 变化 

分 析 图 3 可 见 ,气温 突变 点 前 适 窒 在 东北 三 省 分 布 面积 最 大 的 是 温带 针 闻 混交 林 , 为 55.89 万 km ,主要 
分 布 于 松 嫩 平 原 . 辽 河 平原 .三江 平原 S 小 兴安 岭 和 长 白山 海拔 小 于 800m 左右 的 区 域 ,在 气温 突变 点 后 该 植 
被 类 型 的 适宜 分 布 面积 为 47.22 用 km, 总体 适宜 面积 减少 了 8.67 万 km ,其 中 主要 收缩 区 域 为 辽河 平原 , 面 
积 为 11.71km ,扩张 区 域 为 大 兴安 岭 东 部 和 南部 以 及 小 兴安 岭 海拔 在 400 一 500m 和 长 白山 海拔 在 800 一 
1100m 的 区 域 ,面积 为 3.04 万 km”。 

气温 突变 点 后 适宜 分 布 面积 增加 的 植被 类 型 为 暧 温带 落叶 阔 叶 林 和 温带 草原 ,其 中 暧 温带 落 阔 叶 林 在 气 
温 突变 点 前 的 适宜 分 布 面积 为 55.39 万 km ,主要 分 布 在 除 大 兴安 岭 ,小兴安岭 和 长 白山 海拔 低 于 600m 以 下 
的 其 它 区 域 在 气温 突变 点 后 适宜 该 植被 类 型 分 布 的 面积 为 60.90 万 km ,总 体 适宜 面积 增加 了 5.51 万 km”， 
成 为 东北 地 区 适宜 分 布 最 广泛 的 植被 类 型 ,其 中 主要 收缩 的 区 域 在 东北 南部 沿海 区 域 ,面积 为 1.62 万 km , 扩 
张 的 因 域 在 松 嫩 平原 北部 .小 兴安 岭 海拔 在 300 一 400m 和 长 白山 海拔 在 600 一 800m 之 间 的 区 域 ,面积 为 7.13 
万 RE 

气温 突变 点 后 适宜 分 布 面积 减 小 最 多 的 植被 类 型 为 寒 温 带 针 叶 林 ,也 是 气温 突变 点 前 后 适宜 面积 变化 最 
大 的 植被 类 型 。 在 气温 突变 点 前 适宜 分 布 的 面积 为 46.54 万 km? , 主要 分 布 在 大 兴安 岭 小兴安岭 .长 白山 、 松 
嫩 平 原 北 部 和 三 江平 原 大 部 ,气温 突变 点 后 面积 为 30.91 万 km ,总体 适宜 面积 减少 了 15.63 万 km? ,其 中 主要 
收缩 的 区 域 在 松 嫩 平 原 北 部 和 三 江平 原 大 部 ,面积 为 15.83 万 km2 ,扩张 的 区 域 为 大 兴安 岭 海拔 在 800m 以 上 
的 区 域 ,面积 仅 有 0.20 万 km 。 

气温 突变 点 后 适宜 分 布 面积 增加 最 大 的 植被 类 型 为 温带 草原 ,在 气温 突变 点 前 适宜 分 布 的 面积 为 50.47 
万 km ,主要 分 布 在 松 嫩 平 原 大 部 .辽河 平原 中 部 和 北部 .三江 平原 .长 白山 海拔 低 于 600m 以 下 的 区 域 ,在 气 
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图 例 图 气温 突变 点 后 增加 区 域 ”国王 气温 突变 点 后 减 小 区 域 “轩辕 气温 突变 点 前 后 相同 区 域 


图 3 东北 三 省 气温 突变 点 前 后 植被 类 型 适宜 分 布 区 域 变化 示意 图 


Fig.3 Change schematic of the distribution zone for zonal vegetation before and after temperature mutation points in Northeast China 


温 突 变 点 后 适宜 该 植被 类 型 分 布 的 面积 为 53.46 万 km? ,总体 适 宜 面 积 增加 了 2.99 万 km? ,其 中 收缩 的 区 域 
在 辽河 平原 西部 和 北部 .长 白山 南部 ,面积 为 4.14 万 km ” ,扩张 的 区 域 在 小 兴安 岭 海拔 在 300 一 400m 和 长 白 
山海 拔 在 800 一 1100m 的 区 域 ,面积 为 7.13 万 km?。 

高 山 冻 原 在 东北 地 区 适宜 分 布 面积 最 小 , 主要 分 布 在 大 兴安 岭 海 拔 在 900m 以 上 和 长 白山 海拔 在 1800m 
以 上 的 区 域 ,气温 突变 点 前 的 适宜 分 布 面积 为 104.79km2 , 突变 点 后 减 小 了 27.91km?; 亚 高 山 矮 曲 林 的 分 布 面 
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积 也 较 小 ,主要 分 布 区域 为 大 兴安 岭 海拔 在 700 一 900m ,气温 突变 点 前 的 适宜 分 布 面积 为 3.03 万 km? ,突变 点 
后 面积 减 小 为 2.79 万 km 。 
3.2.4 植被 类 型 分 布地 理 中 心 变 化 


气温 突变 点 前 后 各 植被 类 型 分 布 的 地 理 中 心 均 发 生 了 不 同 程度 的 移动 (图 4) ,从 移动 方向 上 看 ,北部 的 
高 山 冻 原 和 寒 温 带 针 叶 林 向 西北 方向 移动 ,其 它 植被 类 型 分 布 的 地 理 中 心 均 向 东北 方向 移动 。 从 移动 距离 上 
看 ,最 大 的 是 南部 地 区 的 亚 高 山 矮 曲 林 ,移动 了 135.44km; 其 次 是 温带 针 阔 叶 混 交 林 ,移动 了 71.80km; 移 动 距 
离 最 小 的 是 南部 的 高 山 冻 原 , 只 移动 了 1.39km; 其 它 植被 类 型 的 移动 距离 在 22.20 一 71.47km 之 间 。 
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4 气温 突变 点 前 后 植被 类 型 分 布 区 域 中 心 转移 示意 图 


Fig.4 The change of geographic center of the zonal vegetation before and after the temperature mutation points in Northeast China 


3.2.5 ”推算 分 布 区 与 实际 分 布 区 对 比 

由 于 本 研究 缺乏 实地 调查 资料 ,因此 借助 2001 年 和 2013 年 MCD12Q1 数据 的 分 类 结果 ,与 3.2.2 中 各 植 
被 类 型 的 分 布 区 域 进行 对 比 ,在 MCD12Q1 数据 的 分 类 结果 中 落叶 针 叶 林 、 混 交 林 、 落 叶 阔 叶 林 和 草地 分 别 对 
应 3.2,2 申 植 被 类 型 的 寒 温带 针 叶 林 ,温带 针 阔 叶 混 交 林 、 暧 温带 落叶 阔 叶 林 和 温带 草原 。 人 研究 中 首先 利用 
MCD12Q1 数据 的 分 类 结果 提取 2013 年 各 土地 覆盖 类 型 与 2001 年 的 空间 差异 ,增加 部 分 的 数据 分 别 设 为 
A1l1、A2、A3 和 A4, 减 少 部 分 的 数据 分 别 设 为 B1、B2、B3 和 B4; 再 设 3.3.2 研究 结果 中 各 植被 类 型 在 突变 点 后 
增加 部 分 的 数据 为 al 、a2 、a3 和 a4, 减 少 部 分 的 数据 分 别 设 为 bl1、b2、b3 和 b4, 分 别 求 取 Al 和 al、…、A4 和 
ad Bl1 和 bl、…、B4 和 b4 的 共同 区 域 并 计算 面积 。 

结果 表明 (图 5) , 寒 温 带 针 叶 林 无 共同 减少 区 域 ;共同 增加 的 区 域 分 布 在 大 兴安 岭 海 拔 在 800m 以 上 的 地 
区 ,增加 面积 为 0.01 万 km , 占 Al 面积 的 4.82% ,另外 82.57% 的 Al 面积 分 布 在 气温 突变 点 前 后 未 变 的 适宜 
区 域 。 温 带 针 阔 叶 混 交 林 无 共同 减少 区 域 ;共同 增加 的 区 域 在 大 兴安 岭 东部 .小兴安岭 海拔 在 400 一 500m 区 
域 和 长 白山 海拔 在 800 一 1100m 的 区 域 ,增加 面积 为 1.27 万 km?, 占 A2 面积 的 21.46% ,另外 61.01% 的 A2 面 
积分 布 在 气温 突变 点 前 后 未 变 的 适宜 区 域 ; 暖 温带 落叶 间 叶 林 无 共同 减少 区 域 ; 共 同 增加 的 区 域 在 松 嫩 平 原 
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北部 小兴安岭 海拔 在 300 一 400m 和 长 白山 海拔 在 600 一 800m 之 间 的 地 域 ,面积 为 0.78 万 km’, 占 A3 面积 的 
20.23% ,另外 59.94% 的 A3 面积 分 布 在 气温 突变 点 前 后 未 变 的 适宜 区 域 ;温带 草原 共同 增加 的 区 域 在 小 兴安 

岭 海拔 在 300 一 400m 人 山海 拔 在 800 一 1100m 的 区 域 ,面积 为 0.10 万 km , 占 A4 面积 的 8.03% ,另外 
80.15% 的 A4 面积 分 布 在 气温 突变 点 前 后 未 变 的 适宜 区 域 ;共同 减 小 的 区 域 在 辽河 平原 西部 和 北部 以 及 长 白 


山南 部 , 占 B4 面积 的 18.72%, 男 外 70.24% 的 B4 面积 分 布 在 气温 突变 点 后 的 不 适宜 区 域 。 可 见 气 候 变 化 对 
各 植被 类 型 的 分 布 已 经 产生 一 定 影响 


2 寒 温带 针 叶 林 温带 针 疗 叶 混交 林 
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图 5 推算 分 布 区 域 与 实际 分 布 区 对 比 


Fig.S Comparison map between the projected distribution regions andactual distribution regions 


4 结论 与 讨论 


4.1 缚 论 

以 东北 三 省 为 研究 区 ,以 1961 一 2013 年 的 196 个 气象 台 站 逐日 整编 资料 为 基础 数据 ,以 温暖 指数 的 适宜 
范围 为 标准 ,研究 气候 变化 对 东北 地 带 性 植被 适宜 分 布 区 域 的 影响 ,并 结合 2000 年 和 2013 年 土地 类 型 覆盖 
数据 (MCD12Q1) ,对 推 ed 吉 果 表明 : 

(1) 近 50 年 来 大 兴安 岭 气 温 突 变 点 为 1982 年 ,其 它 地 区 为 1988 年 ;东北 三 省 年 降水 量 没 有 明显 的 突 


变 点 。 


(2) 通 过 建立 WW 与 经 度 .纬度 和 海拔 高 度 的 线性 拟 合 方程 ,对 缺少 气温 观测 地 区 进行 数据 插 补 ,再 利用 
差 值 对 持 补 数据 进行 订正 ,订正 后 的 取 值 与 实际 计算 值 更 为 接近 。 
(3) 气 温 突变 点 后 适宜 东北 三 省 生长 的 主要 植被 类 型 种 类 无 变化 ,但 各 植被 类 型 的 适宜 分 布 区 域 均 发 生 
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变化 ,气温 突变 点 前 适宜 在 东北 三 省 分 布 面积 最 大 的 是 温带 针 阔 混交 林 , 气 温 突变 点 后 适宜 分 布 面积 有 所 下 
降 , 其 中 收缩 的 区 域 为 辽河 平原 ,扩大 的 区 域 为 大 兴安 岭 东 部 和 南部 以 及 小 兴安 岭 海 拔 在 400 一 500m 和 长 白 
山海 拔 在 800 一 1100m 的 区 域 ;气温 突变 点 后 适宜 在 东北 三 省 分 布 面 积 最 大 的 是 暧 温带 落 阔 叶 林 , 其 中 的 收 
缩 的 区 域 在 南部 沿海 区 域 ,扩大 的 区 域 在 松 缴 平原 北部 、 小 兴安 岭 海拔 在 300 一 400m 和 长 白山 海拔 在 600 一 
800m 之 间 的 区 域 ;气温 突变 点 后 面积 减 小 最 多 的 是 寒 温带 针 叶 林 ,也 是 气温 突变 点 前 后 适宜 面积 变化 最 大 的 
植被 类 型 ;气温 突变 点 后 面积 增加 最 大 的 是 温带 草原 。 

(4) 与 主要 植被 类 型 实际 分 布 变 化 相对 比 , 在 寒 温带 针 叶 林 .温带 针 阔 混交 林 、 暧 温带 落叶 阔 叶 林 和 温带 
草原 中 ,2013 年 比 2001a 增加 的 区 域 分 别 有 4.82% 、21.46% 、20.23% 和 8.03% 的 面积 落 在 理论 推算 的 增加 区 
域内 ,由 此 可 见 从 理论 推算 和 实际 分 布 上 均 证 明了 气候 变化 对 研究 区 主要 植被 类 型 分 布 可 能 已 经 产生 了 
影响 。 

4.2 讨论 

(1) 本 研究 在 检验 近 50 年 来 气温 突变 点 的 基础 上 ,推算 气候 变化 对 东北 三 省 主要 植被 类 型 适宜 分 布 区 
的 影响 ,并 选取 典型 植被 类 型 利用 现 有 的 遥感 分 类 数据 进行 对 比 , 人 研究 结果 与 书 有 的 研究 结论 有 很 多 一 致 之 
处 ,如 程 肖 侠 等 引 认 为 在 维持 当前 气候 条 件 不 变 的 情况 下 ,东北 森林 树种 组 成 基本 维持 平衡 ,但 气候 增 暧 使 
主要 针叶树 种 比例 下 降 , 阔 叶 树 种 比例 增加 ,与 本 研究 结论 中 气 漫话 变 点 前 后 东部 地 区 适宜 在 东北 生长 的 主 
要 植被 类 型 种 类 无 变化 ,但 寒 温 带 针 叶 林 和 温带 针 阔 混交 林 在 气温 突 受 后 后 适宜 分 布 区 面积 下 降 , 暧 温带 落 
叶 阔 叶 林 适 宜 分 布 区 面积 增加 的 结论 基本 一 致 ;在 本 研究 中 气温 突变 点 后 温带 针 阔 叶 混 交 林 的 适宜 分 布 区 域 
已 经 侵 和 人 大兴安岭 东 部 ,这 一 结论 验证 了 程 肖 侠 等 :5 在 研究 赤 兴 安 岭 森林 演 替 中 指出 的 气候 变化 将 可 能 使 
大 兴安 岭 会 以 温带 针 阅 叶 混 交 林 为 主 的 预测 ; 冷 文 芳 ' 引 在 HADCM2SUL 和 CGCMI1 模型 预测 假设 100a 后 全 
球 平均 气温 和 降水 量 增加 的 基础 上 ,认为 寒 温 带 针 叶 林 主要 建 群 种 (兴安 落叶 松 和 云 冷 衫 ) 的 覆盖 率 有 明显 
下 降 , 这 一 结论 与 本 研究 结果 一 致 ,但 对 于 本 研究 中 的 温带 针 阔 混交 林 的 主要 树种 红 松 ,本 研究 的 结论 更 倾向 
于 CGCM1 方案 的 预测 结论 ; 曼 寒 冰 ':55 等 认为 气候 变化 使 兴 安 落 叶 松 和 和 白桦 的 分 布 向 西北 漂移 ,这 与 本 研究 
中 寒 温 带 针 叶 林 适 宜 分 布 区 域 的 地 理 中 心 向 西北 移动 的 结论 一 致 。 

(2) 本 研究 得 出 的 气候 变化 对 研究 区 主要 植被 类 型 分 布 可 能 已 经 产生 了 影响 的 结论 ,是 在 理论 推算 和 遥 
感 分 类 数据 相 结合 的 基础 上 得 出 的 这 当 申 主要 误差 来 源 有 .一 是 理论 推算 方面 ,本 人 研究 推算 气候 变化 对 植被 
类 型 分 布 的 影响 是 建立 在 气候 变化 突变 点 基础 上 的 ,气候 变化 突变 检验 的 结果 是 近 50 年 来 具有 温度 因子 发 
生 了 突变 ,因此 认为 温度 因子 在 项 中 起 了 主要 作用 ,尽管 在 本 研究 中 已 检测 到 降水 无 显著 变化 ,但 降水 .自然 
灾害 、 人 为 活动 等 诸多 因素 对 精 被 类 型 分 布 也 具有 一 定 的 影响 ,因此 本 研究 得 出 的 结论 具有 一 定 的 片面 性 ,但 
方 精 云 . ”和 分 春 静 等 '* 认 为 , 单 因 素 标准 对 于 确定 大 区 域 植被 带 划 分 具有 和 较 好 的 指示 意义 ,而 且 刘 向 培 "3 
的 研究 结果 表明 待 东 北 地 区 温度 对 地 表 植 被 动态 变化 的 影响 大 于 降水 ,因此 本 研究 只 根据 温度 因素 推算 气 
候 变 化 对 植被 类 型 分 布 的 结果 仍 具 有 一 定 的 参考 意义 ;二 是 遥感 分 类 数据 的 精度 问题 ,本 研究 直接 采用 了 美 
国 NASAALPDAACEOS 数据 中 心 提 供 的 土地 类 型 覆盖 数据 (MCD12Q1) ,由 于 缺乏 历史 地 面 调查 数据 ,因此 无 
法 对 划分 类 精度 进行 检验 ,分 类 误差 也 会 影响 本 研究 结果 准确 程度 ,虽然 在 讨论 中 中 已 有 相关 研究 成 果 支 持 
了 条 研究 的 结论 ,但 本 研究 给 出 的 每 种 植被 类 型 分 布 区 域 变化 的 具体 数据 还 有 待 进一步 验证 。 

《937 虽然 物种 分 布 与 气候 因素 有 密切 关系 ,但 物种 的 迁移 对 于 气候 变化 的 影响 往往 有 一 定 的 滞后 性 ' ， 
也 会 受到 自身 的 迁移 能 力 人 为 干预 和 诸多 外 界 条 件 的 影响 ,而 本 研究 在 推算 植被 类 型 分 布 适宜 区 域 变 化 时 
并 未 考虑 上 述 因 素 的 影响 ,这 也 是 把 推算 结果 与 植被 类 型 实际 分 布 变化 作对 比 时 ,共同 增加 (减少 ) 面积 占 实 
际 增加 (减少 ) 面积 比例 较 低 的 原因 之 一 。 
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